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ABSTRACT 


locally similar solutions for boundary layer on a 
circular cone in ass3nrimetric conical flow is analysed in 
spherical polar coordinates. A series solution in the 
ascending powers #f angle of attack a and thickness para-r 
meter ^ , that is the ratio of local boundary layer thick- 
ness to the l*cal radius of transverse curvature, is 
developed. The various terms of this series are split into 
universal functions, which are independent of prescribed 
conditions namely, semi-vertex cone angle and the free 
stream conditions. The leading term is the well known 
Blaussius Equation and the remaining equations are linear 
and form a recursive set of equations. These equations 
are integrated numerically using Runge-Hutta Method, Though 
only first order terms in a and are considered in this 

rV" 

analysis, it can be extended to consider higher order terms 
in a and Tc , 



LIST OE SYMBOLS 


Chapman- RuLisen constant 
coefficient of friction 

stream function related to meridinal velocity, 
defined by Equation (2.1.26) 

stream function related to meridinal velocity, 

defined by Equation (2.1.21) 

stream function related to circumferential 

velocity, defined by Equation (2.1.21) 

universal functions related to meridinal velocity 

universal functions related to circumferential 

velocity 

total enthalpy, non-dimensional, defined by 
Equation (2.1.6) 

differential operator, defined by (2.3.20) 
differential operator, defined by (2.3.28) 

Mach Number 
pressure 

distance from vertex of cone 
local transverse radius of cone 
non-dimensional distance from vertex of cone 
Reynold number 
local Reynold number 
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T Temperature 

u meridinal velocity 

V velocity normal to the surface of cone 

w circumferential velocity 

X perturbation in meridinal velocity due to angle of 

attack in inviscid flow over cone 
y perturbation in circumferential velocity due to 

angle of attack in inviscid flow over cone 
z perturbation in density due to angle of attack in 

inviscid flow over cone 

SUBSCRIPTS MD SUPERSCRIPTS 

0 conditions at surface of cone at zero angle of 

attack in inviscid flow 

e conditions at the edge of boundary layer 

- denotes dimensional variables 

^ differentiation with respect to similarity 

variable 

&REEK SIMBOBS 
a angle of attack 

p defined by Equation (2,1 .15) 

y ratio of specific heats 

6 ■ boundary layer thickness 

6* displacement thickness -as so dated with u 



displacement thickness-associated with w 
similarity variable, defined by Equation (2.1.21) 
semi-vertex cone angle 

angle between radius vector and cone axis 
defined by Equation (2.1*6) 

perturbation in pressure due to angle of attack 

in inviscid flow over a cone 

coefficient of viscosity 

thickness parameter, defined by (2.1.20) 

density 

shear stress 

azimuthal angle 

stream functions, defined by Equation (2.1.13) 



CHiOTER 1 
IHO?RODUGTION 

Motivation and Scope of the Problem 

When a slender cone is placed in supersonic stream, 
the inviscid flow over the cone is conical in nature, that 
is flow properties do not change along a ray from the vertex 
of the cone. This holds for the case of a cone at an angle 
of attack also. If the fluid is viscous there is a boundary 
layer formation on the surface of the cone. The conical ihvis- 
cid flow affects the detailed structure of the flow in the 
boundary layer and provides the necessary boundary conditions, 
at the edge of the boundary layer. At an angle of attack, 
pressure on the wind-ward side is higher than the pressure 
on the leeward side, and thus there exists a circumferential 
pressure gradient. Since near the surface of the cone, the 
inertia of the fluid is less than that on the edge, fluid near 
the surface tend to move in the direction of pressure gradient 
more closely than the fluid near the edge. Hence low energy 
fluid near the surface is drained from windward side to leeward 
side. Thus thickness of the boundary layer is more on leeward 
side than on the windward side. In the forward region of the 
cone, the thickness of the boxmdary layer is comparable with 
the transverse curvature of the cone. This is more so on f | 
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leeward side of the cone. Thus considering the locally 
similar boundary layer, it is possible to predict the 
results more accurately near the vertex of the cone and 
one can go far down the leeward side. 

We know that boundary layer grows thicker as the 
Mach number in the external stream increases. At hypersonic 
speeds, boundary layer remains no longer thin on most par+ '' 
the cone in comparison to the transverse^ curvature of the 
cone. Hence, thickness effects should be considered in any 
analytical development of the boundary layer flows, at high 
Mach numbers. 

From the above physical consideration, it has been 
found necessary to consider the locally similar boundary 
layer and estimate more accurately the detailed structure 
and the overall characteristics of the boundary layer flow 
over the cone at angle of attack. 

T’rom mathematical point of view, we hope to gain a"^ 
advantage over the existing approaches to solve this problem. 
Most of the boundary layer analysis is carried out in the 
cylinderical polar coordinate systan, while for the cone 
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problem, spherical polar coordinate system seems to be more 
natural. Because, all the inviscid flow analysis are in sper- 
ical polar coordinate system and this coordinate system elimi- 
nates the need of using Mangler's transfoimation, to convert 
the equations of motion to almost Cartesian foim. In the present 
investigation, boundary layer flow is analysed in sperical 
polar coordinates. A similarity transformation is found for 
this coordinate system and a series solution in the ascending 
powers of the angle of attack, a, and thickness parameter r. 


I 
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Cone angle 

distance from vertex of cone 


u> developed Jlppendix 1 shows that use of spherical polar coordi- 
nate system yield same result for thin boundary layer. The 
various terms of this series are split into universal functions, 
which are independent of prescribed conditions namely, cone 
angle, angle of attack, free stream condition etc. The leading 
term is nonlinear and remaining equations form a recursive 
set of equations, which are integrated from "ttie stirface to 



4 


tlie edge of the boundary layer by 'using Runge-Kutta method 
on I.B.M. 7044 Computer with extended precision facility. 

The values of these functions are tabulated upto six decimal 
places without rounding off the values. The effect of the 
prescribed conditions appears in the coefficients of these 
universal functions. Thus once these tables are available, 
computation of various dynamic and thefmo dynamic characteris- 
tics of the flow can be done by simple algebraic manipulations 
These universal functions are tabulated in Appendix 5* 

The computed results. for two particular cases show 
that as thickness parameter increases, there is an increase 
in skin friction along the generator. But there is a decrease 
in contribution of incidence in skin friction as thickness 
parameter increases. Expressions for coefficient of friction 
and displacement thickness are derived and can be used for 
the particular case. There are certain limitations of this 
work also. These limitations are as follows: 

(i) Viscous-inviscid interaction has been neglected, 

(ii) Though the boundary layer is considered thick, 

.0-momentum equation has been neglected. It does not 

affect the analysis upto the terms of order a, but 

2 

once terms of order a or higher are considered, 0- 
. momentum equations should also be considered. This 
equation is 

0p ' ■ 2 

Got 0 
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(iii) Results upto terms of a & Y only are e-valuated, 

that is only two terms of the series have been consi- 
dered. Within this approximation, the pressure gradien'^.- 
is favourable in the circumferential direction. As such 
one should not expect to approach the separation lime 
with the present analysis. In order to include the 
terms of order ^ and higher, Taylor expansion of Sin 0 
should have terms of similar order in (0 - 

(iv) At large angle of attack, near the zone of separation, 
vortex is shed into the inviscid zone, thus modifying 
the inviscid flow around the cone. In this case this 
analysis will not hold. To date, no analytical or 
experimental work is available to predict the inviscid 
flow, taking into account the shedding of the vortex 
into the inviscid zone. 

Review of the Past Work 

In reference (1) Maccoll has fitted isentropic coni- 
cal flow to a eonical shock, cone being at zero incidence. 

Due to conical nature of the flow and circumferential symmetry, 
the three dimensional equations of motion are reduced to one 
of ordinary non-linear boundary value problem. This equation 
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is to be solved numerically. Boundary conditions are provided 
by Rankine-Hugnoit relation across the shock, ^t the shock 
and on the surface normal velocity is zero. 

In reference (2) Kopal has tabulated, the numerical 
solution to the equation given by Maccol (1), for various 
semi-vertex cone angle and various free stream Mach manber* 

In this table, the values of radial and normal velocity com- 
ponents are presented as function of angle 6 for a given 
value of the semi-vertex cone angle, 0 , and Mach number, 

, These tables also give surface velocity on the cone, 
shock wave angle and drag coefficient as a function of cone 
angle and free stream Mach niunber. 

In reference (3) Stone has dealt with the problem 
of cone at small angles of attack in supersonic stream. The 

effect of the angle of attack is considered to be a pertur- 
bation over the basic solution of Maccoll (l) for the case of 
zero angle of attack. In this reference the linear terms in 
angle of attack a are considered. The resulting equatioiB are 
simplified using the equations of zero incidence case. These 
resulting equations are to be integrated numerically. The 
uniqueness of the solution to the approximation taken is also 
proved, under the assumption that flow properties and their 
derivatives depend continuously on the angle of attack a. 
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The equations are integrated starting from the conditions at 

the shock wave for the symmetric flow until the conditions 

on the cone are satisfied, that is normal velocity component 

■becomes zero for 0 = 9^. 

c 

In reference (4) Kopal has tabulated the numerical 
solution of the equations given by Stone (3). The coefficients 
of a in; velocity components, u, v and w, pressure and density 
are tabulated for a series of semi-vertex cone angle and free 
stream Mach numbers. It has been found that axial force compo- 
nent does not change from its value in symmetric flow but 
there exists some nomnal force component. This is zero in case 
of zero incidence. These axial and normal force coefficients 
are also tabulated in this reference. 

In reference (5) Stone has extended the analysis of 

reference (3') to include the terms of second order in angle of 

attack a and derived the equations governing the coefficients 
2 

of a in the expansions of velocity components u, v and w, 

pressure and density. But in deriving these relations, terms 

3 

of the order of a and ~ are neglected. Since normal velocity 
V becomes zero near the surface, the equations are not valid 
near the surface of the cone. Boundary conditions are also 
derived in similar manner. These equations contain singularity 
in the sense that few terms have normal velocity v in the 
denominator, which goes to zero as cone surface is reached. 
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This singularity has been taken care of, by making certain 
substitutions. 

The numerical solution of equations derived by 
Stone '(5). are tabulated by Kopal in reference (6), In these 
tables various coefficients are tabulated as a function of 6 
for various values of the free stream Mach number and semi- 
vertex cone angle. The coefficients in expressions of velocity 
are refered to maximum velocity while coefficients of pressure 
and density are refered to their corresponding values at zero 
angle of attack. These tables also give corresponding shock 
wave angle, axial force coefficient and normal force coefficient. 

In reference (7) Perri has improved upon Stone’s (3) 
result upto terms of order a and has pointed out that the form 
of the solution chosen by Stone does not satisfy entropy con- 
servation equation at the surface of the cone. Since cone sur- 
face forms a stream surface, the entropy on the surface of 
the cone is constant but Stone’s solution gives variable entropy 
on the surface of the cone. In this reference asymmetric conical I 
flow has been studied and it has been shown that there must ' 

exist singular points in the plane of symmetry. At these singu— j 

lar points entropy is discontinuous. This singularity exist | 

at the point where the circumferential velocity has negative ; 

circumferential gradient, Perri has shown that in most of the | 
regions stream lines are radial and near the surface of the 
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cone these stream lines move circumferentially and converge 
at the leeward plane of symmetry. Since each stream-line has 
different entropy, there is singularity in entropy at leeward 
plane of symmetry. Also there is a region near the surface 
of the cone, where the entropy varies singnificantly. This 
region has been termed as ’Tortical’ layer by Perri, This 
vortical layer affects the derivatives of velocity component 
of zero and first order near the surface of the cone. The term 

, which were taken to be zero in Stone’s work, are of the 
order of a in the vortica,! layer. In this layer it affects 
the value of , by quantities of the order of a and pressure 
by the quantities of the order of 4 - , It has also been 
shown in this reference that for determination of the pressure, 
complete Bemaulli equation should be used. The pressure dis- 
tribution calculated in this manner, compares well with exper- 
imental analysis. 

In reference (S) Young and Siska developed the equ- 
ations for flow properties in a region bounded by a cone at 
a large angle of attack and its attached shock wave for the 
use of tabulation of Ecpal (2,4»6) , .In Reference (4,6), the 
coefficients of a in the expansion of the velocity components, 
pressure and density are tabulated, but the boundary conditions 
are formulated for the conical flow field at zero angle of 
attack case. Thus while using Kopal table this point need be 
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talcen care of* The equation to correct the values of flow 
properties are derived in two ways, one hy assuming that each 
intermediate conical surface still remains conical and is 
rotated by angle of attack oc, and the other by assuming 
that intermediate conical surface are distorted linearly in 
a shape from that of solid cone to that of shock wave at an 
angle of attack. These equations are derived with respect to 
the wind axes. Equations for the cone surface characteristics 
e.g. 0^ etc, are also derived with respect to the body 

axes. Thus this reference improves upon the v^rork of Kopal (4j6), 
so as to satisfy boundary conditions more accurately. 

In Zopal tables (4,6), results are tabulated with 
respect to wind axes, and Young and Siska (8) have given formu- 
lae to use those tabulations in a proper manner the axes still 
being the wind axes. Only cone surface characteristics are 
derived in body axis system. In reference (9) Robert and ^ley ; 
have given the formulae similar to the one given by Young and 
Siska, but with respect to the body axes. Thus using Zopal table ; 
one can get flow properties at a point whose coordinates with : 

f 

respect to the body axes are known, Eor the derivation it is \ 

'U j 

assumed that the intermediate surfaces are yawed with solid | 

^ l 

cone and remain circular cones. The relations for first order I 

[ 

and second order terms in a are separately derived. Computed [ 

results show good agreement with experimental result even for 5 
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a > e . This has been pointed out that there is an error in 
c 

tabulation in the values of the perturbation in circumferential 
velocities. 

^correct ~ table) - (2 x/SinS) 

In reference (10), Willett extended the analysis of 
Perri (7) to the terms of order of a . Perri (?) gave correction 
to the solution given by Stone (3), Kopal (4) for the presence 
of vortical layer. Willett has given the equations for correction 
to Stone (5), Kopal (6) results, due to the presence of vortical 
layer. In this reference components of velocity on the surface 
of the cone u and w, pressure and density are evaluated, using 
Kopal (4,6) and Robert and Riley (9). Thus the flow properties 
are obtained in body axis system. These relations are derived 
by assuming that entropy on the cone surface is constant and 
equals to the entropy of the fliiid in windward plane of symme- 
try. It is also assumed that pressure distribution is correctly 
given by Stone (5) to the second order. In the derivation of 
these relations it is assumed that normal velocity is zero. 

Hence formiiLae are valid on the surface of cone only and not 
in the flow field. Experimental results show that the entix^py 
variation with angle of attack is correct upto a » 20° for 
semi-vertex cone angle 0^ = 15° and Mach number 3,53* These 
results on the surface of cone can be used as the boundary 
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conditions at tlie edge of 'boundary layer on the surface of 
the cone. 

In reference (11) Bulakh has extended Stone’s (5j5) 
results inside the vortical layer. The solution inside vor- 
tical layer upto the accuracy of a is o'btained. This solu- 
tion reduces to the Stone’s results outside the vortical 
layer. It is also shown that derivative&of entropy, radia,! 
and circumferential velocity with respect to polar angle 
0 becomes infinity at the surface of the cone. In 

Stone’s (3,5) solution there is a logrithmic singularity 

2 

in entropy if series is terminated upto a , but if series 
is carried upto infinity i.e, all the terms in a^, (n- a nat- 
ral number) are considered, such logrithmic singularity 
disappears. This is shown by taking general solution of the 
entropy conservation equation near the surface of the cone. 
Thou^ logrithmic singularity is shown to disappear when 
all the terms in a are taken, any numerical work will still 
have this problem. This reference also proves the assumption 
of Willett about the pressure on the cone and the entropy 
Jump anal 3 rtically* 

In reference (12) Cheng solved the problem of 
hypersonic flow around the cone at an angle of attack. 
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CJheng considered two parameters, one the limiting density 

K-1 

ratio across shock wave » and other being the angle of 

attack. The form of solution is same as that of Stone (3,5) 
except for terms corresponding to the aa©* parameter density 
ratio, . Cheng gave analytical solution upto the first 

order terns in density ratio and angle of attack. Expressions 
for circumferential velocity and pressure are given upto 
second order approximation also. Glie derivative of the lift 
coefficient with respect to angle of attack agrees well with 
Kopal ( 4 ) values. Analytical expression for entropy distri~ 
button near the surface of cone is given, first order correc- 
tion to pressure, density and radial velocity are given by 
assuming that the singularity does not affect leading appro- 
ximations of normal and circumferential velocity v^w. These 
corrected results show that the effect in pressure is in 
third order only and in the normal and circumferential velocity 
it is of the second order. 

In reference (13) Munson solved the problem of the 
cone at an angle of attack using the technique of matched 
assymptotic expansions. Outer expansion is equivalent to \ 
Stones (5 ) expansion. The inner expansion, which is valid 
in thin vortical layer near the body describes flow proper- 
ties in the vortical layer* In the inner expansion -coordinate 
0 has been changed to 0“ in order to take care of logrithmic 
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perturbation in density. These two expansions are matched 
assymptotically and then combined to give uniformally valid 
solution throughout the flow field bounded by cone and attac- 
hed shock wave. This uniformly valid;, solution is given both 
for first order and second, order terms separately. It is 
also shown that pressure does not change inside the vortical 
layer. He also shows that upto the second order the entropy 
singularity still remains on the surface of cone and does not 
go into the fluid. 

In reference (14) Babenko et.al discussed numerical 
method of solving the Euler's equation with proper boundary 
conditions, A difference scheme for solving the problem of 
three dimensional flow has been set up. The stability of the 
method of solving these difference equations is also consi- 
dered, The method is applied to conical flow equations and 
equations are solved numerically for various values of angle 
of attack, Mach number, and semi vertex cone angle, Eor solv- 
ing the equations for a given angle of attack, first the equa- 
tions are solved for zero angle of attack, and then given a 
small angle of attack Aa, computation is repeated by taking 
zero angle of attack solution as basic data. This process is 
repeated till the required angle of attack is achieved. This 
process also required some movement along the generator to 
establish the conical nature of the flow. The results are 
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tabulated as a function of ^ = (6-9^)/(0-9^) , where 9g shock 
wave inclination and the azimuthal angle, for various Mach 
numbers, angle of attack and semi vertex angle. Inclination 
of shock wave for each azimuthal angle is also tabulated. 
Calculation also shoV(?s that when Mach number and angle of 
attack are kept constant and semi -vert ex cone angle is inc- 
reased, shock wave is nearest to the cone in the windward 
plane of symmetry and later becomes nearest to the cone in 
the leeward plane of symmetry. The results of computation 
show that for small angles of attack variation of velocity 
component with respect to ^ is' nearly linear but as angle 
of attack is increased it becomes more and more nonlinear, 
while in the analysis of Stone (3,5), this dependence is 
same for all angles -of attack. These tables can be directly 
used for finding the flow variables for a given tabulated 
case. But for untabulated case one has to interpolate the 
results in order to get the flow properties and it is diffi- 
cult to say what interpolation rule should be used. 

In reference (15) Moore has derived the three dime- 
nsional boundary layer equations on cone for conical flow. 
These equations were derived using boundary layer equations 
in cylinderical polar coordinates and then transforming,-'. these 
equations using Howarth •transformation, Mangier 's transfo|-- ; 
mat ion and Blaus siuja transformation. 
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In reference (16) Moore has analysed the laminar 
boundary layer on a right circular cone at an angle of attack 
in a supersonic stream. The outer solution to inviscid flow 
is taken from Kopal (4). Moore considered that for small 
angles of attack, the boundary layer quantities differ# slig- 
htly from the corresponding quantities in the flow in case 
of zero angle of attack. These perturbation quanties are 
assumed to be of the same -form as the outer inviscid flow. 

The system of equations governing these perturbation quantities 
is solved to yield the effect of small angles of attack on 
the velocity profiles, skin friction and boundary layer thick- 
ness etc. Actually the rate of change of these quantities 
with respect to angle of attack are evaluated at zero angle 
of attack. It has been shown that if we take the case of 
a/0 » 1 , we get the equations for the compressible viscous 

flow about a yawed infinite cylinder. 

In reference (17,18) Moore has analysed the problem 
of laminar boundary layer on a ri^t circular cone at large 
angle of attack to a supersonic stream, in the plane of symme- 
try only. It is shown that in the windward plane of symmetry 
boundary layer becomes thinner when angle of attack increases, 
but in the leeward plane of symmetry analysis fails to give 
unique solutions except for the small angle of attack. Beyond 
a certain critical angle of attack boundary layer flow does 
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not exists all. The ratio of this critical angle of attack 
to the cone serai-vertex angle is given as a function of Mach 
number. This non-existence of solution on the leeward side 
indicate a possible separation in the leeward plane of symme- 
tiy. 

In reference (19) Probestien and Elliot considered 

the problem of transverse curvature effects in compressible 

axially symmetric laminar bpundaiy layer flow. In this reference 

the thickness of the boundary, layer is assumed to be compara- 
radius of 

ble to the^transverse curvature of the body. The transverse 
curvature effects are characterised by the parameter which 
is essentially the ratio of the local boundary layer thickness 
to the local transverse body radius. Modified Mangier ’s trans- 
formation is used to change the boundary layer equation with 

terms 

transverse curvature ^o almost Cartesian form. The analysis • 
is mainly concerned with the case of transverse curvature 
parameter being nearly or less than unity. The velocity ^d 
temperature are expressed as a power series in transverse 
curvature and problem is solved for the case of power law 
bodies, both compressible and incon^iressible flow and exponen- 
tially growing body. The zeroth approximation gives the Mang- 
ier ’s result and higher order approximation modifies the 
result to include the effect of transverse curvature. The first 
order analysis for the case of cone and cylinder with zero 
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pressure gradient ah.OY'js appreciable change in skin friction 
and heat transfer rate. This change keeps on increasing v\?ith 
Mach number at the edge of the boundary lawyer. 

In reference (20) Reshotko extended the analysis 

of Moore (17) for the case of non- insulated surfaces. The 

boundary layer equations are especialised to the plane of 

symmetry and are solved numerically. In the plane of S 3 nanetry 

equations become ordinary differential equations. In this 

analysis boundary layer is assumed to be thin compared to tra- 
radius of 

nsverse^curvature of the cone and ra.tio of specific heats is 
taken to be constant. While solving the differentia,! equations 
in leeward plane of symmetry, it was noted that except for 
very small angles of attack either the solution is indeter- 
minate or non-existant , It was found that for the case of 
Prandtl number Pr = 1 , velocity and temperature profiles 

(ifi — fTT ) are same. The circumferential velocity profile exceeds 
^ o'-^w 

the external circumferential velocity. It is noticed that 
heat transfer increases significantly with angle of attack 
for a given free stream Mach number and cone semi-vertex angle. 

Por free stream Mach number of 3.1 and cone semi -vertex angle 
012=5° » heat transfer is more than twice at, an angle of attack 
8° as compared to heat transfer at zero angle of attack. The 
heat transfer coefficient changes faster for almost insulated 



19 


case than for cooled waP. case. In thia reference method is 
given to extend the ana,lysis to very large angle of attack 
using the results of case of yawed infinite cylinder. 

In reference (21 ) Vvedenskaya solved the prohlen of 
infinite circular cone at an an^e of attack, in a viscous 
heat conducting fluid. Using finite difference method, self- 
similar solutions of boimdary layer equations, are computed. 
Boundary layer is considered thin compared with ^transverse 
curvature of the body. In thro reference after similarity 
transformation, simila.rity variable is scaled so asto have the 
value unity at the edge. The equations are first solved in the 
plane of syimnetiy and then in the meridian planes away from 
the plane of symmetry. Equations axe considered to be initial 
value problem in the circumferential variable. The values of 
the variable at the edge are taken from Babenko et , al (14), 

It was found that for small angles of attack solution of equa- 
tions was smooth over the whole surface and the solution matched 
with that solution in the plane of symmetiy in the leeward side. 
The boundary layer also grew monatonically from windward plane 
of symmetry to leeward plane of symmetry. Heat transfer and 
friction decreased moixotonically with (j). But for larger angles 
of attack (a >2,5° for *= 10° 5 a > 5° for 9^ = 20°) smooth 
solution was not possible and circumferential derivatives of 
velocity components u,w and temperature T did not go to zero 



at leeward plane of ssnunietiy. Thus there was a discontinuity 

in first derivative at leeward plane of symmetiy. In such cases 

it vi/as not possible to solve equations in the leeward plane 

of syrnsetry. The displacenent thickness on leeward plane of 

syrmetry is several times larger than that on the windward 

plane of syinmetry. Por = 5, and temperature of the body be-- 
ok 

ing 500 j displacement thickness changed from 0*8 on wind- 
ward side to 5-*1 on leeward side, for 6^^ = 10° and a = 7*5°. 

Tor the same cone when angle of attack was 2,.5*^.it changed 
from 1.2 to 3.3* '^he circumferential derivative of heat tran- 
sfer first increases with a and then decreases at leeward plane 
of ssnmmetry# 

In. reference (22) iwyer computed compressible laminar 
boundary layer flow over a sharp cone at an angle of attack 
slightly greater than scAi— vert css ~o one angle. The boundary 
layer equation after using similarity^ Howarth transformation 
are solved niimerically. These equations are solved as an initial 
value problem in circumferential variable, and boundary value 
problem in slmilcxlty variable Boundary layer is consi- 

dered thin compared to the^transverse curvature. Tor boundary 
conditions at the eclge, the result of Munson's (13) work has 
been taken, i'or initial condition in boundary layer equations 
were solved in windward plane of symmetry by quasi-linearicciion 
technique. Then boundary layer equations were solved by the 
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finite difference technique developed by Blottner (24), after 
some modifications. Nonlinear teims were evaluated at previous 
station in <}). Momentum and energy equations were decoupled 
between (j) grid locations byevaluating dependent variables 
at the previous station. The continuity equation was written 
in integral form and solved numerically. Steps in ({) were redu- 
ced till convergence in velocity and enthalpy was obtained. 
Computed results show that in most of the regions relative 
heat transfer that is actual to zero angle of attack case 
agree for small a with experimental value if second order expa- 
nsions for inviscid flow are used. But as the angle of attack 
is increased the experimental value deviates from computed 
value. It has been found that beyond a certain angle of attack 
first order inviscid theory gives serious errors. It is also 
shown that as angle of attack increases beyond the cone semi- 
vertex angle, second order expansion of Munson (13) gives 
inadequate boundary conditions, as the boundary layer structure 
comes out to be incorrect both ^qualitative and quantitative way 
on the leeward side. 

In reference (23) Boerick has done work almost similar 
to Dwyer (22), In reference (22), reference quantities were 
arbitrary. While in this case radial and circumferential velo- 
cities are normalised with respect to their respect ive ■ edge 
velocities. Enthalpy is also normalised with respect to the 
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edge value. In reference (22) circumferential coordinate is 
the actual angle <|), while in this case it is (j) Sin 0^. In 
the leeward plane of symmetry, in order to avoid the numeri- 
cal difficulties due to vortical singularity, the circumfer- 
ential coordinate is changed so as to eliminate the singularity 
from the equations. In the plane of symmetry equations were 
solved hy the ’method of accelerated successive replacements’ 
proposed hy Lieherstien (25) and applied to boundary layer 
problems by Lew (26) and Storm, (27)* Blottner method was used 
for solving equations away from plane of symmetry. The compu- 
ted results were within 6% of the experimental value for rela- 
tive heat transfer, when calculations are based on the experi- 
mental pressure distributions. In the leeward plane of symm- 
etry error is more than that at other meridinal planes. 





CHAPTER 2 
&EIfERAL ANALYSIS 


2,1 Eormulat:' . ■ n 

Boundaiy layer equations in spherical polar coordinate 
are as follows: 


^ ( r^ Sin e p u ) + 1^ Cr Sin e p v) + |» (r p w ) * 0 
dr ^ 

( 2 , 1 , 1 ) 

- 8u V 3u _w gu _ ^ 1 gp, 

^ gr~ r r Sin"? 5^ r p dr 


P 


1 

r Sin e 


3 

30 


sin 9 5|-) 


(2,1,2) 


u 



V 8w 


w 8w . 

r sin e r 


1 

P r Sin 0 



+ 


_1 3 

P r Sin 0 80 



Sin,0 



(2.1,3) 


«/^Jl t0M 

and for the case of Prandtl No. Pr = energy equation hecomes 


H, 


1 

2 



+ 



L 


(2,1.4) 
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Boundary conditions 

at u=0 = t= w 

at edge u u » w -♦ w 

© 0 


( 2 . 1 . 5 ) 


I’on-dimensionalising the variables with respect to the 
inviscid flov^ quantities on the surface of cone at zero 
incidence. 



and using linear viscosity law 




0 


0 



( 2 . 1 . 6 ) 


( 2 . 1 , 7 ) 


where 

Uq - velocity at the surface of cone at zero incidence -in 
inviscid flow. 

n - density at the surface of cone at zero incidence in 
inviscid flow. 
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- temperature at tlie surface of cone at zero incidence 

in inviscid flow* 

P UqI 

Re-r = — 5 Reynold number 

Mo 

1 - representative length. 

Gq - Chapman Rubisen constant 

the equations (2.1.1) to (2,1.4) and boundary conditions 
(2,1.5) transform to 


(R^ Sin ePu) + 1^- (R Sin ep v) + (Rpw) 0 

( 2 . 1 . 8 ) 

n ^ -L I. ^ ^ = » 1 9-E. 

^ aR R 3© R Sin 9 ^ R P 9R 


_2 9 Q 

2 (ip Sin 9 qq) 

pR^ Sin 9 ao 


(2.1.9) 


aw . V 9w . w 
" 5 80 R Sin e 


aw ^ ^ ^ 

^ R p R Sin 9 ^ 


+ 


1 

P R^ Sin 9 


a_ 

a© 


(|- Sin 9 
o 


aw \ 


( 2 , 1 . 10 ) 


H. 


2 2 
u + w + 


2V 


E 

P 


(2.1.11) 
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Boundary conditions 

at 0=0 u=0=T = 'w 

( 2 . 1 . 12 ) 

at edge u -♦ u j w w 

© 0 


Using stream functions f and X. such that 

Ijr = Sin e p u 
0 

Xq = Rp w 
and 

= - R Sin 9 v 


(2.1,13) 


and taking Taylor’s expansion of Sin 9 around 9^ upto two terms 
Sin 9 = Sin 9^ (1 + Cot 9^ (8~9^)) 

(2.1.14) 

= Sin 9^ (1 + p 0) 
c 

fd . 

where p = Cot (2,1*15'^ 

Equations (2,1,8) to (2,1,11) and boundary conditions (2, 1,12) 


transform to 
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fo 




fR Sin JPR Sin 9^ 


PR Sin 


V® (1-P ®) f x@ 

+ — 

C-. /-» • . 


?Q- (1-P 0 

:2nrr 


]piF^n 9^ U pR^ Sin 9^ /pR" Sin 9^ I ^j) 


4^ 




= ■ p IS 

4' 


T 


, , 1 1^0 

(l+p0) ^ 2 


/R^ Sin e~ \ e 


(2.1.16) 


1^0 (X© , _ ,X®, , _X 6 ,Xe 


rt ^ 


pR** Sin 9 q ^ pR'^ Sin 9^ 


^ - -.K — <7 rAlSLl J- - 

rr— ^(jRh * si„ 0 ^PR^I) 


V® x® 


v2 ^4 


1 9 

3 


p^- R‘^ Sin 9^ PR Sin 9^ 


1 3 \ T 




pR^ 


-^j^r- (1+p0) (^)o| (2.1.17) 


and 


H. 




0 


7R^ Sin 9, 


(1-pS)^ + +1^ ^ (2.1.181 
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Boundaxy Oondltions 


at 0=0 


at edge 


^0 “^0 ^ 

I 

’■>- 7 ^ (1-P0) V. 

0R‘^ Sin 0^ 


Xe 

^R 


( 2 , 1 . 19 ) 


w 


e 


It is known that the boundary layer -thickness varies 
as the reciprocal of the square-root of Reynolds number and 
in case of cone transverse curvature varies directly with R, 
the distance from the vertex taking the projection of both 
the boundary layer thickness and the radius of transverse 
curvature on the plane normal to the axis of cone, one has 


boundary layer thickness 
radius of transverse curvature 

Qq fsL. 

Jr v 



( 2 . 1 . 20 ) 


Row making similarity transformation from (R, 0, <})) to 
(^^>7?, <t» such that 




1 R C^) 

2 p. 


o 


^ /pd 0 

0 






and ^ 

f = Sin 0^ f 





( 2 . 1 . 21 ) 


One gets from the Equations (2.1,16) - (2,1 .18), after making 
use of conical nature of outer inviscid flow, the following 
equat ions 


- - 2 . g ({) hn _ _ j ^^%I2Lht±. - 2 2 

^ ^ 7 } ~ 3 Sin 0^ 5 Sin 0 p ^ Sin 0 3 



[ils? 


23/2 ^ 

+ ^ 3 ) . mr w 


V sS-8^ ^ fp;' 
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^ i d7?+ £i (E^) 

2 Sin V o p ^ 2 Sin 9^ 


1/2 


O 3/2 (P.'/'/ 1 c 2 )'^!fZLk>-(?i) (E~) 

^3'* Sin 9^ So'’ i p ^ 2 Sin 6 ^ S ^ 


1/2 




2.V2_2,p,Va j;, |(E_)Va_£2, 


-(f)' 4(^) 


0 /> 


3 "P, 


. /| f, 4). 


-r) 'p^jf 


0 


(2.1.22) 


- f g 


2 ^(|) ^7)7? 1_ 


7777 . 3 Sin 9 ^ 5 Sin 9^ 


^ 'm 


2 J2I 


^ 3 Sin 9 ^77 ^ 


2 - 2 ^ L ^ L %?nl. ^ I f. g 

■^3®T)V3 PSlne^ S) 5 3 i 




1/2 


g 


77 ^ 


(^) 

^0 


1/2 g7777 

p 


7777 


0 (2,1.23) 


(H„ - f?, - gt - J:^ ?r) + 


2 22^ P^ 

■7? - s?|"Y~i p 



0 O 

I'i^ 


A 


\ 


(2.1.24) 
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using (2. 1 ,21 ), -■boundary conditions (2.1.19) ‘becomes 
Boundary conditions 

at 7? = 0 f = f = = g({) = g =, 0 

at7>-“ Jjj (1--^ j- ^ {E_) / 


(2.1.25) 


In Eq.uation (2,1,22) and Equation (2,1.24) terms 
2 

of the order ^ are neglected 'because ^faylor series expansion 
of Sin 9 is taken upto two torns, that Is 


Sin 0 = Sin 9^ (1 + Cot 9 (9 9„)) 

C Q c 

= Sin 9- (1 + p 0) I and from ('2.1.21) 
c 



0 invqlves the terms in ^ and hence, unless the hi^er 

order terms in Taylor series expansion of Sin 9 are considered, 

\ 

one- should not consider the hi^er order terms in ^ in the 
above equations. 

Since the "boundary condition on f as Tjl ^ is a compli- 
cated one and difficult to satisfy while solving the equation 
numerically, the following substitution is used in order to 
simplify the boundary condition. 
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1/2 


V 




(2.1.26) 


substituting (2,1,26) into Equations (2,1,22) to (2,1,24) 

% 

and neglecting terms of ^ , one gets 


. 2 6(1) . g . ^2 &nH di _ 2 ^2 

y sin e^, Tpf “ 3 Sin 6^^T}* 3 Sin 


V2.n,V2, 1 

f__ /t 

o P 


7? 




P 


•<r- 




o 


1 


/2 !i (E_)^^^ 7 ii - gl / 

J3 Sin 9^ p -' p 3 Po 0 


1/2 

TT 

:= 0 


(2.1.27) 


2 S(t . 1 g ^(D „ 2 ^7) S7?j/ 

“ ^ ^7?" 5 SlnT^ ^7?” J Sin 9^ P ^7)7?*^ T Sin 0^ 


2 ^ 4> ^ 2 ^d) 

j St^XTj + 5 pSin 0 


g 

c7)7}7? 



7? 

g (f / 


dl^ 

“ 


7? 


- { - f"7) 


V dTj 
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= 0 

2,1.28) 

“ ^ 7 } ® 7 ) y " p (2,1.29) 

Using (2.1,26) 'boimdaxy conditions become 

at 7} = 0 f ^ g = f = f . = g . = g = 0 

' ^ (2.1.30) 

at 7 ) -♦ 00 fjj-* xLq 

S77 



2.2 Series Expansion of Stream Functions 

Assuming that flow variables f , g and P can be expressed 
as a power series in ^ , the thickness parameter, one can write 

P= Po +^Pi + i%2 

/ ^2 ( 2 , 2 . 1 ) 

S = §0 +^§1 + §2 

Herep^, f^, g^j p^, f^, g^ etc. are functions ofTj and 
alone. 

Substituting (2.2,1) in (2<,1.,27) to (2,1.29) and separating 
the coefficients of like powers of ^ one gets following equation^. 
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^0 ^07/77"^ ^ ^07777"^ 3 Sin 0 ^ p 


% ^oriv ^ 


2 ^ \ 2 „ ‘2- , ^ ^ 
3 Sin e" ^ ^07} (j) T ^077+ ^0777777" 


r 


^o ^^Tfr}~) ^3 


\ 9 V2 ^ 1/2 '' . 

-.(t) 


(f) Cp^) ^ 0777)0 " ^ ^°-nv^ ^ 


,2^5/2 .r, , 1/2 


( 2 . 2 . 2 ) 


“ 3 Sin 0 ^ ^^ 1 ( 1 ) ^07777*^ ^17777 ^o(|)^ ~ 3 Sin 0 ^P ^o 


717) 


^‘'777)^ 3 Sin 0 Q ^^177^077^ ^ ^077 ^■‘ 7 ? 

'- 2 m 


2 (E_)^^^ f /n dTj- /T-" 




'T?777?o 

V2 2\n ,V2 p 1 


Po 


■17;7?7) 




2 ^ -- 1 

^0 ^07777' 3 Sin 0Q ^o(|) 3 Sin 0^ p ^0 ^07777 


3 Sin 0^ ^077^0 77 ({> “ T ^ o7 )-^ o7 } 


SrtTj ■* 


2 

T 


1 


P„ Sin e,, 


■" «o7,JJ7?= ° 


( 2 . 2 . 4 ) 
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P 3/2 ^ 1/2 


) o 


L 


’irp} 


■° *7775 0 P< 


,2,3/2 .r, ,1/2 


7 ? f 




£_ / s 1 ^d) / 

3 Sin e^ ^®1(j) ^oypfj'^ ®o4 ^ 17 ^ 7 ^ ^ 5Sin 9^ -p ^17)7} 


^7)7} 

+ IT ? i dr)+ g _ ) - I II, 


■'■ ^1 3 Sin 0^ ^®n7^^l7J(j) ®177^077<1>^ 3* 


V 




1 


3 ^P 


o 


°V °7) 0 /3, 


‘■07) "17? ^ , 3 siST,, ^ 


■ ®1 7)7)7" ? ^07) ®1?7 ~ I 3 '■t 

/r,E_,V 2 go 7)77 _ 

/ 5 P„ 




0 


H 1 _ ?JUr ^ ^ 0 

®0 ^07)- SoT) p 


"O 


^ 077^177 ■^® 07 ;^ 177 “ Y ~1 


y 


= 0 


02.2.5) 

(2,2.6) 

(2.2.7) 


Boundaiy conditions 

at 7) = 0 fjj = f ,| = = f,,^ = gg = g,| = goi() " 


0 


” ^OT) " "e > ^ 17 ? 

«^ 7 ) " ’ 817 , - 0 


81 $ - « 07 ,' 8 i 7 ?= ° 

( 2 . 2 . 8 ) 
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Now Equations (2.2.2), (2.2.4) and (2.2,6) sho^d be 
solved for the case of thickness parameter ^ being zero. In 
order to include the effect of ^ , Equations (2,2,3), (2,2.5) 
and (2.2,7) should be solved usihg the solution of Equations 

(2.2.2) , (2.2.4) and (2.2.6). Equations C2.2,6) and (2.2^7) 
are to be used to eliminate and from Equations (2.2.2), 

(2.2.3) , (2.2.4) and (2.2.5). 

Outer inviscid flow for the case of small angle of 
attack is given by 

u = 1 + ax Oos (j) 
w = ay Sin ^ 

® (2.2.9) 

5= 1 + az Cos (|) 

^ = 1 + aACos (j) 

Expanding the stream functions and density f^, g^, f ^ , g^ , 
and in the similar form one has 

^0 = ^oo + ^ f^^ 

g^ = ay Sin (() g^^ (2.2.10) 

Po = Poo + ^^01 


I 
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and + OCX Cos ([) 

= ay Sin (() 

^1 ' Ao ^11 


( 2 . 2 . 11 ) 


Here ^oV^ol’ ^oo’ Pol’ ^10’ ^11’ ^11’ 10^^11 frictions 

of 7} alone. Prom Equation (2,2.9) one gets 


E— = 1 + aXCos ({) 

Pq 

H(h 

and ^ = - a A Sin <|) 

2 

neglecting the terms of oc and higher order. 


( 2 . 2 . 12 ) 


Substituting Equations (2,2. 10) and (2.2.12) in Equations 

2 

(2.2.2), (2.2.4) and (2.2.6) after neglecting temns of a 
and higher order and collecting the terms of like powers of a 
one gets 


f f + f =0 

^o© oojijj oorjTjT) 

^ 007 ^ 77^10 


(2.2.13) 

2 y g f 

J X Sin 0^ ^o 1 007777 

(2.2.14) 


g, 


’o1 


+ f 


vm 


00 





2 .xf 

J y:i Sin 9^ 



(2.2.15) 
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H 


o 



gV £o 

y-'' Poo 


0 


and 

y Pq 

Y*’’ P 


X 

X 


z 

X 


+ fo1 0 

Poo 


(2.2.16) 


(2.2.17) 


Boundary Conditions using (2,40) in (2.38) gives 


at = 0 

f = f =0 

oo 007^ 

at 7 ] - °o 

"^007? ^ 

( 2 . 2 . 18 ) 

atl) ^ 0 

0 

li 

0 

H 

0 

at ?7 


( 2 , 2 . 19 ) 

V 

at 7 } ~ 0 

Sol =.Sol 7 J= ° 

at 77 — 


(2,,2..20) 


K'ow substituting Equations (2,2.10), (2,2,11), (2,2,12) and 
(2.2,13) into (2,2.3), (2.2.5) and (2.2.7) and collecting 
the coefficients of like powers of a, one gets 


%yiv ^ 


^oo'fj ^10 


rj 


^ 2 „ 
+ y £ 


00 


m 


10 


1 


^ /V 

'00 i . n 

rfr)7f ^oo 


? 5/2 

d?7 + (y) f 


7? f 
/ 


Poo 



(2.2.21) 
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o1 iOTi ^3^ ) ^4x 00 , 
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• 'I 5 V_;' 3-'>* <f>''l^'-w '•JitC*’’ 
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+ { Poo ^ °\rf ^ P , 


'n i . 0.3/2 . ?foo P 01 


777)1 

<1^ 

^ d7? 


7m 


00 


- ^ f f - 


5 o1^ 10 3 Sin e^ x 


fe-s. -i *^*00 S-11 + fio„o,Sol) - 3 ^107,^01 
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vv 


'71 ^'7) 


2 1 ^ ^oo7)7} 



z ^oOT^ol 
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T )2x p 


00 


A 


00 


Poo 
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( 2 , 2 . 22 ) 


®''^7}' ^^°°77«117? 


3/2 


= - (^) f.n Sc 




^/x ■ 
a( 7? 

5\ ' f o 

/ m - f f.n g. 


'r -00 " 01 ^^ Poo ' 5 10 -01 


r ) 


3/2 

■" ®al7)^“®7){poo 


P, 


^7} '■ 2 A Pq ':iO . 2 _ „ 

^ + 3 ^ y Sin e^ ^ 3 ®o1^ lo^j 


- r S . 


V ^ _ /r Jaizp?^ 


(2.2.23) 


'V Po P 10 

___ ^ 


■00?^ 01^ 7-1 p 


= 0 


00 


(2.2.24) 
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y 

+ ^01^^10^ y-1 

„ ^10 

00 

r-i + 2 Eu. 

^ ^Poo 

^ ^o1 ) 
^ 0 

i-'OQ 




(2.2.25) 

Boundary Conditions 




TJcing end 'Z 2 .S) 

one gets 



at 77 = 0 f^Q = f^Q^= 0 

at?7'» 

0 

0 

(2,2.26) 

atT9=0 f^^=f^^^=0 

at?7 ^ 

fj,^ "*■ 0 

(2.2.27) 

at r) = 0 0 

atT)"*’ 

^11 Tj"* ^ 

(2.2.28) 


These equations are ordinary differential equation of 
third order. These are to he solved numerically in order to 
find the solution to original problem, 

2.3 Reduction of Stream Ttinctions into Universal Functions 

Since the equations governing the stream functions 
^o1' g^1> ^-{0* ^11’ ^11 non-homogencous and non-homogenious 
part contains the terms like x, z, A etc., which depend upon the 
particular flow situation, these need he solved for each parti- 
cular case of free stream mach number, semi-vertex cone angle 
and angle of attack. In order to avoid this problem of solving 
these differential equations each time separately the stream 
functions are broken into universal functions. The differential 
equations governing these universal functions are independent 
of the particular flow situations, and their linear combination 
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gives these stream functions. In this linear combination the 
coefficients of the universal functions depend upon the parti- 
cular flow situations. Thus the problem is to be solved only 
once and any particular case can be derived by talcing the proper 
linear combination of these universal functions. 


Lett ing 


'OO 


= ]? 


oo 


(2.3.1) 


from 


(2,2,13) and (2,2,18) one gets 


puf + p p« - 0 
00 00 00 


( 2 . 3 . 2 ) 


Boundary conditions 


atV ^ 0 B = I" = 0 

' 00 00 

at 7? - « I” -* 1 

' 00 


(2.3.3) 


(Here prime denotes differentiation v;ith respect to ), 
Brom the definition of p one has 


P = 




At the edge of boundary layer for zero incidence one has 


Po = ^ 

P U 

00 0 


1 


(2.3.4) 


y m2 
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and from (2,2.16) one has 


Ho = 1 + 


<1- 

(Y-D^mF 


(2.3.5) 


where is the Mach nttmber at the edge of the boundary 
layer in the case of zero angle of attack. 


Now letting 


y-i m2' 


Sc1 = ^OO + (1 G 


o" ^011 


(2.3*6) 


from (2,2.15), (2,2,16) and (2.3.6) one gets 


ttt If off 

'o11 + ^00 ®o11 ■ T ^00 ®011 


2 

3 


00-^ 

(2.3.7) 


Boundary conditions 

at 7) » 0 G, 

at 7) “ G 


o11 

t 

o11 


^011 = 0 


0 


(2.3.8) 


letting 

f 


“ T X Sin e„ ^00 + 3 " X Sin 9)^011 


Z. 


o1 


-22. 


3x Sin 8, 


(1 + M^) (2.3.9) 
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and after substituting in eq.uations (2.2.14) and (2,3.6) 
one gets 


ttt 

^o11 

^00 ^o11 ^00 ^o11 ~ ° 


Boundary 

conditions 


at 

7] = 0 ®o11 

yoii = 

at 

7) - “ ®'o11 ■* 

1 


(2.3.10) 


(2.3.11) 


and 


^o'ia ^00 K-\2 * ^00 ^o12 


- C ®o11 


(2.3.12) 


Boundary conditions 
at 7? - 0 

at 7) ^ 

Similarly, letting 
^10 ~ ^101 


^o12 "" ^o12 


5’ — 0 

^o12 


0 


+ (Y- 1) ^102 


(2.3.13) 


(2.3.14) 


and su-bstitutlng in (2.2.21) and using (2.2.16) and (2.?. 2) 
one gets 

+ ^00 ^" o 1 + y ^00 ^ f C ^ 1 o 1 


.2 3/2 « 


ri 


i rs.'n - §■ (2.3.15) 
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2 5/2 

'0 ^oo '"oo 


» „ o 5/2 « V 

( K . + ®n 1 l) + C|) ^00 { ^00 


u 


2.^/^ aT, - 1 (C - 


+ ^o11,' o 00 


o1V 1o1 


+ (|)^^ 7) Poc, (Poo ®o11^ 1^101 ^^00 + ®0l4 




Boundaiy .,^'onditions 
at 7) = 0 
at 7) -* =» 


Gj J. -1 •=» ^ ^ 


G- 


^111 

I 

111 


111 

0 


(2.3.22) 


Ii (^112^ ” 


- I (|)" ' *'oo ^^00 + ®0l 1 


3/2 


00 


tl ^ \ 

) (7) - 1 S’ -B) 

^ ; \fi ■^Qo ^Qo 00 


/2s^/^ B g” + ~ ■‘‘ ^o11^ 

- (3) ”2 ^00 ^o11 ^ 2 ^ 00 oil 


( / ^00 '^^^00 ^00 ^o11 { ^00 


2? 


. 2 5/2 » 

1/2%' „ r 1 


o O 


- 1 p<r,o (p„„ + ®oii^ ■ y ^ioi ®oii 


3 102 '00 


1 ,22/2,' (p' * 0 ' ,) (7?- P„o <0 - C) 


- f ®oo ^"00 ” '' 0 II 
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+ 1 C|)% <„1 + (F< + <ii) 


-) r t i po~. »f ttf . , f \ 

3 ^1o1 ^o11 " 2J 3 ^00 ^o11^ “ ^oo ^00 “ ^oo 


^/j <11 - j/Kii 


1 [2 p " 

2 j3 oil 


(2.3.23) 


at 7) ^ 0 


at 7] - 


Cti12 - ~ ^ 




(2.3.24) 


I (S1I3) = - I Poo <11 (T?- ^00 <0 - « 


+ 5- (f) ^on 0 - Pno) <117- rPmO ^ol 


0II Q 00 


5 "1o2 "oil 


+ I <| 6 ' ■ P;o <11 ( V - Pio Poo - ^00) + §■ ^00 ^012 


3/2 , 


2 ’ ‘ 

^ Tn m 


+ 5 <11 <02 - iff <11 ^00 <0 - <o) 


1 12 n 


00 00 00' 


tJj <11 


(2.5,25) 
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Boimdary conditions 

at ?7 = 0 
at 7] -» oo 


&113 - - 0 


^113 ° 


(2.3.26) 


^11 “ ^111 X Sin 0 


(Y- l)Ma 


4 


+ ■'’ “o ^114 * 3-sii-e; ^"5 


X 


A 


^116+ ""o ^117 


(2.3.27) 


and 


tff 


n 1 * - , 2 -ra ” 

(f ) = f + I’oo ^ + T ^00 ^ 3 ^00 


(2.5.28) 


and 


using (2.3.27), (2.2.22), (2.3.6), (2.3.9), <2.3.14), (2.3.19), 
(2.2.16), (2.2.17) etc. one gets, 

3/2 


(^111) = - ^o11 <01- <11 ^ <0 ^o1l5 


75 /p T) 

+ (|)^^^ <0 ^ ^011 ^v - 1 ^ioi <11 ■ y <°i 
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Boundary conditions 


at?)=0 = = 0 


at 7] 


^111 ’ ° 


( 2 . 3 . 30 ) 


ft 


p u p 5/2 

l2 (P112) = - I - I-oo) - (t> V ®oo <^011 


" _ 1 _ ' 

■’ 0 I 2 "^ 00 ' 


+ ^..-10 "• ^00) ~ (3’) 7 } ^00 ^^o1t '*’ ^o 12 '" ^00^ 


2,^/2 ^ m «t m . 

■ 7 } (^o11 " ^00 + ^012^ 


2,5/2 ^ M (t 


tf 


V 


+ (f) - Poo) f ^00 


2 5/2, Tl 

+ (|) P^o, { '^011 + ^012 - ^ 00 ) -ill- 


■ ? ^1o1 ^^o11 ■'■ ^o12 " ^ 00 ) " y ^00 ‘'HI 


y ^00 ^ 1 o 1 “ y ^ 1 o 1 ®o 11 ■ 9 ^ 1 o 1 ^o 12 " ^oq) 


^ 3/2 ir tl tl . 

- (y) ^^o11 ^o12 "" ^ 00 ^ 


( 2 . 3 , 31 ) 
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Boundary conditions 
at T) = 0 

at Tj-* °° 


^’112 " ^112 


^112-* ° 


(2.3.32) 


1 ,2n^/^ ^ r-n- b b'. - bI^) 


(Pl13> = - 5 ”o2 *'o11 ■ ? ^“■'1 


■00^ 


^ S'n11 <-V- ^00 ®O0 - ^00^ 


" 3jJ 00 -^011 


ft! 


^|T^o11 ^00 ^00 


00 - ^ 00 ^ jjl ^00 { ^00 ^on'^n 


nt V 1 » 


1 

+ ^ 




V 


2 T?" r B (1 - 5'^^) dT) 
f ^o11 i "^00 ^ 00^ I 


7) 


ilfj'o / 0 - *oo) 

3/2 .// 


+ 5 I 3 ■‘^00 Q ■^o1 1 

V 


t t 


- (|)'"X'o" ■'" ^00 ^oo ^011 I ^011 ^102 


»l 


i B B - i T “ 

“" 3 ^o11 ^^102 2 3 o11 


i "Z"" < t If 

"-/Tv ^OlV ^00 

(2.3.33) 


Bb undary c ondit ions 

at 77 = 0 

at 7? “ 


I'll? = ^113 “ ° 


B ■* 0 

^113 


(2.3.34) 
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~ 


|<o2 ^^011 ^012- ^00^- T ^012^101 


1 /2 


3/2 


j ( 3 ) 77 ^o 12 


’00 


- i ^00 ®00 ®00 - ^oo> ^*'011 


tl t1 V 

- P + 


00 


o12‘ 


3 X 3 . 


. 9 3/2 t» 

j (j) 7) ^00 ^o12 

^ ( 7 )- ^00 - (^011 - ^00 ^012^ 


-VI§vK,2 


. (X- ^00 - (/"1 ' ^ 0 ° 


h'oo </l 1 - /o - /12) 


7) 


+ 


1 ^o12 “ ^00^ { ^00 " ^00^ 


+ 3 (j) 


3/2 , 


0 I 2 0 


V 

f 1 d7] 

00 ' 


7) 


1 F** / S’ 'i'n 

+ 5 - ( 3 ) ^00 j 012 
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1 (L\ -p / (1 

• 2 : v ’ s ;/ -^nn ' 


3 "3’ 


^l.b (SnU + ^012 - ^00) 


_ s" I P„o Too ^^011 ■ ^00 


'00 ^ 00 


~ 1^101 <12 - 1 C ^112 -5 ^011 ^101 


- I ^1o2 ^^o11 ^o12 ■ ^00^ 


2 ^1. T?”_ + G_.. J'lo ^ !lo J 


“ ^ ^1o1 ^o12 


- f (^00 ^102 ^o11 ^1o2 + 


f 9 ^/2 H 

I '^011^^012-^00)+ "o 12 


- ("I - ^o^ '^^11 + ^o12 - ^00 ) 


P 3/2 „ I t - p' ) (2.3.35) 

+ (I) J'oo ^00 <^0l1 + ®012 ^OO’’ 


Boundary conditions 

at 7? = 0 


at 77 


^114 =* ^114 " ° 


^114 ° 


(2.3.36) 


1 . 1 T. KANPUR 

CENTPAt tJBRARY 
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"■ 


^012 ^'lo2 " 9 ji ®00 ^012 (V- ?oo - 


p" ) 

00"^ 


? 5/2 , 

(^) S' jP 
^ 00 o12 


^00 I’oo - C) 


1 TT? *” ^ r If 

-5/? *'o12 (V- Poo ®oo - Foo> 


1 ri" ? I 

rjT^oo /I'oo ^012 17? 


1 2 V2 I. ^ ^ 

+ 5- Cj) P„12 / »oo (1 - Po|) a?? 

1 p 3/2 n ^ 

+ 5 <3 ^ ^00 { *’o1 2 (-^ - /f ) 17? 

C {®oo /o 177-1^012^102 


? /o <>113 - ?/l2 /o2 - 111^012 C1 - Poo) 


1 2 t t 

+ i >* F P P 

5 3 00 00 o12 


(2.3.57) 


Boundary conditions 


at 77 = 0 

at 77 ■*■ ®° 


” ^4+C — 0 


115 ’ 115 


/l:5 “ 0 


(2.3.38) 
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^1 ^ z jj ^00 -5-^00 ■[ ^00 7 l§ ^00 ( 2 . 3 . 39 ) 


Boundary conditions 


at 7J — 0 — 0 


at 7? 


116 ~ "^116 


^116 “ ° 


(2.3.40) 


and 


I, (F117) = + ^ (f) ?oo ^co ^00 - *' 00 ^ 


+ ?,/f C {\o (1 - ^00) 17?+ f /f 'C (1 - 


(2.3.41) 


Boundary conditions 


at 7^—0 — 0 


•■"117 “ ^117 


(2.3.42) 


at 7? — 


F. 


117 


These equations (2.3.1), (2.3.6), (2.3.10), (2.3.12), (2.3.15), 
(2.3.17), (2.3.21), (2.3.23), (2.3,25), (2.3.29), (2.3,31), 
(2.3.33), (2.3.35), (2.3.37), (2.3.39) and (2,3,41) are to be 
solved recursively with proper boundary conditions. Once these 
equations are solved with proper boundary conditions, one can 



get f^Q 5 g^^ and and lienee f and g for a 

given particular flow situation, thus giving detailed behaviour 
of the flow. 

These equations are solved numerically by Runge- 
Kutta Method for solving third order differential equations 
as given by Collatz (28). 
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GHAPTER 3 
RESULTS 


The various houndary layer quantities e.g. skin friction, 
houndary layer thickness etc. can be evaluated by using 
Equations ( 2 . 2 ^ 1 ), ( 2 . 2 . 10 ), ( 2 . 2 . 11 ), ( 2 . 3 . 1 ), ( 2 . 3 . 6 ) , ( 2 . 3 . 9 ), 
( 2 * 3 . 14 ), ( 2 * 3 . 19 ) and ( 2 . 3 . 37 ) and the universal function 
tabulated in Appendix 3 . The coefficients of these universal 
function are to he taken from Kbpal ( 2 , 4 ). Figure 3 to 12 
present the quantities related to velocity profile for two 
cases, one for semi-vertex cone angle of 5° and free stream 
Mach number of 8,4925 and other for semi-vertex cone angle 
of 7 . 5 ° and free stream Mach number of 8 , 0589 . 

Velocity Profile 

Velocity profile for the case zero angle of attack are 
shown in Figures 3 and 4 for the two cases. Both these figures 
show that as thickness parameter ^ increases, velocity pro- 
file becomes more and more steep near the surface of the cone. 
This increase in steepness implies increase in skin friction. 
This increase in the skin friction is also depicted in 
Figtires 7 and 8, in which shear profile is plotted for the 
two cases. This increase in the shear stress is due tjo decre- 

t 

ase in local Re3niold*s number as increases, thus increasing 
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viscous forces in comparison with, inertia forces, 

Yelocity profile in circumferential direction are 
plotted in Figures 5 and 6 for the two cases. Both these 

graphs shov7 that circumferential velocity profile becomes 
less steep with increase in thickness parameter ^ « This 
implies that skin friction in circumferential direction 
decreases with increase in ^ , This decrease in the circum- 
ferential shear stress is also depicted in Figures 9 and 10, 
These figures show the shear profile of circumferential 
direction. 

The rate of increment in the meridinal velocity profile 
due to angle of attack is given by 



This quantity is plotted in Figures 11 and 12 for the 
two cases. These graphs show that as thickness parameter 
increases, there is decrease in this rate of increment. 

Skin Friotion - 

The circumferential and the meridinal components of 
the viscous shear stress at the cone surface are given by 
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Of^ = r— §-2- (/^-|f0->l0=e^ 
ofr = - 3 — 

Po "o 

Making ike quantities non-dimensional on tlie ri^t hand side 
of these equations one has 



^ 2 Re^ 
r 

1 


= 

^0 

g^,|( 0 ) 


2/^ A 

</ . 2 Re^ 

1 


Cf, = 

J^o 

f-|<| Co) 

Where 

t:)* 

0 

0 



Re^ = 

^0 




Using expansions of f and g and making use of tabulated 
functions one gets 


- 2^ Ri; 


Sin $ |l.0836 + 0.5070 (V-1) 


- 0.8430 C - 0.2049 ()^~1 ) ” 0.0283 
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Cf^ = 2 




C 2 

^ £ 0*4696 ^ 0*4989 0*2301 (K— 1 ) M^s, 


A A 


+ ax Cos (() ^0.7046 4 - 0*1930 ^ ^ 0.0130 ^ 

.+ 0.1150 (V'-D 1^^0.2487 e_ 


+ 0.0460 ^ Sin" 9 “o + 0.5042 € 

o 

-°-9596 - ^ 111 e; • - O-ISSS “o ? 

- 0-2«i rjE&e: “o 


- 0.0254 X si — 


Both of these expressions reduce to result equivalQit 
to that of Moore (16) for ' = 0 and that of Prohestien and 
Elliot (19) for a = 0 and to that of Mangier for oc 


Boimdary lever Ihickness 

The houndary layer thickness, when expressed as 
displacement thickness, associated with the welocities w 
and u are as follows: 
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>(}) 


,0 - .^)ae 

o n 


. J4)ae 

° Pe'^e 


Expressing the quantities in the non-dimensional.- form on the 
right hand side of these equations one has 

"1/2 oo ^ u 

«) f Or. -=^w)a/j 


(j) ~J3 i 




and 


2 °o 


1 


^r ~ JJ Re 


r ^0 


(£ / A - 2^ 

» 0 

o 


Using the expansions of f and g that is for u and w and 
integrating with the use of tables one gets 


- F ° 

“is Ei- 


|o.6518 - 6.87154 + (0*5707 -2.2929 4) 


- 0.52374 (X-1) 


/2 j:o_ 

'3 Re^ 


£ 1.2168 + 1.9218 (*^1) BJq + 



6& 


j- -0.5618 + 0.3754 (f-l) + 0.3367 ] 

+ axCos (}) 1^0.6214 + 0.4924 (^-1) 

- 0.5031 (V-1) ~' x"sin ' e (^-1)^ 

c 

- 0.6376 ^(i?-1) - 4.9899 gf| -'- Q Ci'-I ) 

c 

- 0.0251 (/-1)^ 4.1295 ^rsll“e 

0 

-°-1243 M® + 0.5+80 ^ 

+ 0.2636 ^ (y-1) - 0,6941 ^ gin 9 . ' 0*5505? 

- 2.7985 F- Slii^e " * 5.3916 |t0.9609 (i^-1) 4" 

c 

- 0.1683 (-K-l)^ j] 

These expressions reduces to the results equivalent to that 
of Moore for ^ = 0, 

It should he noted here that this analysis is valid 
only in the limit of vanishing angle of attack a and thicknes‘=’ 
parameter ^ • ■Actually the terms multiplied hy (x and v. 
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represent tlie- rate of changes Vi/ith respect to these 
quantities, evaluated at a = 0 and^ » o respectively. Whether 
or not the absolute change for a small finite angle of 
attack a and thickness parameter^ can be obtained from this 
theory depends on the relative magnitude of the effects so 
computed, rather than on the size of angle of attack a and 
thickness parameter and depends further on the second and 
higher drder terns in the expressions of f , g etc. 



REFERENCES 


l!accoll, J.W,; ^The Gonical Shock Wave Formed by a Cone 
Moving at Speeds’, Proc. Roy* Soc. A. 159 

(1937) pp 459. 

Kopal, Z. ; ’Fable of Supersonic Flow Around Cones’, 

M.I.T,, Pept. of Electrical jEngg, Center of Analysis, 
Technical Report 1, Cambridge, Mass. 1947. 

Stone, A,H«; ’On Supersonic Flow Past Slightly Yawing 

Cones I’, Journal of Mathematics and Physics, Vol*27, 
1948, pp 67-81. 

Kopal, Z. ; ’Table of Supersonic Flow Around Yawing Cones’, 
M.I.T,, Dept, of Electrical Engg, Centre of Analysis, 
Technical Report 2, Cambridge,' Mass. 

Stone, A.H, I ’On Supersonic Flow Past Sli^tly Yawing 
Cones II’, Journal of Maths, and Physics, Vol«30, 
1951, pp 203-213. 

Kopal, Z. I ’Table of Supersonic Flow Around Yawing Cones’, 
M.I.T,, Dept, of Electrical Engg,, Centre of Analysis 
Technical Report 3, Cambridge, Mass, 

Ferri, A,; ’Supersonic Flow Around Cones at an Angle of 
Attack’, N.A.C.A. T.R. 1045, 1951. 



63 


8 . Young, G.B.W. and Siska, G.P,; ’Supersonic Plow Around 

Cones at Large Yaw', Journal of Aeronautical 
Sciences,' Vol, 19, 1952 , pp. 111 . 

9. Roberts, R.C, and Riley, J.L, 5 'A Guide to the Use of 

M,!,!. Cone fables'. Journal of Aeronautical 
Sciences, Vol. 21, 1954, pp. 336, 

10. Willett, J.E. j 'Supersonic Plow at the Surface of 

Circular Cone at an Angle of Attack’, Journal 
of Aerospace Sciences, Vol. 27, I960, pp. 907* 

11. Bulakh, B.M . 5 'Supersonic Plow Around an Inclined 

Circular Cone’, P.M.M. , Vol. 26, 1962 (Journal 
of App. Maths. Mechanics, pp . 460), 

12. Cheng, H.K. | 'Hypersonic Plow Past Yawed Circular Cone 

and Other Pointed Bodies', Journal of Pluid 
Mechanics, Vol, 12, 1964, pp, 169» 

13 . Munson, A.G. ; 'The Vortical Layer on an Inclined Cone’, 

Journal of Pluid Mechanics, Vol. 20, 1964, pp.625* 

14 . Babenko, K.I.; Voskresehkiy, G.P,; Lyubimov, A.N. and 
Rusinov, V.V,; 'Three Dimensional Plow of Ideal' Gases 

Around Smooth Bodies ’ , Israel Program for 
Scientific Translations, Jerusalem, 1968* 



65 


22. Dwyer, H,A, | ‘Boundary Layer on a Hypersonic Sharp 

Gone at Small Angle of Attack’ , A, I. A. A, Journal, 
Vol. 9, Beb. 1971. 

23. Boericke, R.R, ; ‘The Laminar Boundary Layer on a Cone 

at Incidence in Supersonic Blow’, A.I.A.A, Journal, 
Vol. 9, March 1971. 

24. Blottner,' P. ; ‘Finite Difference Methods of Solution of 

Boundary Layer Equations’, AIAA Journal, Vol, 8, 

Peb. 1970. 

25. Lieberstien, H. | 'Course in Numerical Analysis’, Harper 

and Row, Hew York, 1968, 

26. Lew, H.; ’The Use of the Method of Accelerated Successive 

Replacements for the Solution of Boundary Layer 
Equations', A.I.A.A, Journal, Vol, 6, May 1966. 

27. Storm, C.R, ; ‘Application of the Method of Nonlinear 

Simultaneous Displacements to General Three- 
Dimensional Stagnation Point Boundary Layer 
Equations’, A.I.A.A. Paper 68-786, Los-Angles, 
Calif., 1968. 

28. Collatz, L. ; ’The Numerical Treatment of Differential 

Equations’, Springer-Verlag, Berlin, Heidelberg, 
New York, 1966, 



66 


APPEIfPIX 1 

JUSTIEIOAPION EOR USIBG SPHIRICAL POLAR COORDINATES 

Spherical polar coordinate system is used in place 
of the usual cylindrical coordinate system in the present work. 
In order to justify the use of spherical polar coordinate 
system and corresponding similarity variable analysis of thin 
boundary layer on a cone at zero angle of attack is performed 
in coordinate systems and coefficient of friction is compared. 

Boundary layer equations in cylinderical polar co- 
ordinate system are as follows 

(A.1 .1) 
(A. 1.2) 
(A.1, 3) 

It is assumed that Prandtl number is unity and no heat transfer 
takes place. Outer flow is assumed to be conical in nature. 

Boundary conditions 

at y=0 u = T= 0 

at y -* 6 ^ 


^ ^) + (Pr u) = 0 

0X ^ 

^ ^ ^ A— Ch^) 

02 dy p 9y 9y 

■■ ii = Const = H 


(a. 1.4) 
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Non-dimensionalising the variables 


R = 



H = 



X = x/L, y = y/L 


Re., 


'o 



(A.1*5) 


where 1 is some given length parameter and is Chapman- 
Rub is en constant, one gets 


I 3 P ( /OR u) + R V ) = 0 


u 


du 

0X 




9u 

9y 



H = 


(A. 1.6) 


(A.1.7) 

(A.1,8) 


Boundary conditions at 

y =0 u=v =0 

y -*■ 6 u -♦ 1 • 

Using stream function. |r such that 

fy = p. R U 
= -p R V 


CA.1.9) 


(A. 1.10) 



68 


satisfying equation (A. 1.6) identically, one gets 

- (■§) (■§)y = ^ I ^ |y 


(a. 1.11) 


Boundary conditions 
at y = 0 

at j -* 6 


= %= 0 


ry TX 




(A.1.12) 


Using combined Howarth., Mangier and similarily transformation 


X p 

;c = / dx 

^ o 


7? 


J 

f R pdy 
0 


Y2f 




and 


1|J = ^Tiy f(?7) 


(A. 1.14) 


one gets 


f f + f =0 

rm VVl 


(A.1.15) 


Boundary conditions 

atT7 = 0 f = f =0 


(A.1.16) 


at 71 


1 
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Now 


T = 


0^ = 


9y 


T 


0. 


2 e e 

Pe ^ 

,f'3 

i X ; rst 

J f Jj 

'f r 

fyj-TjCo) 

5; — f 

„Co) ■ 


2 ^ 8u 

— 


ly = 0 


(A.1.17) 

(A.1.18) 


(A,1 .19) 


Boundary layer equations over a cone in sperical polar 
coordinate system are as follows* 


9 


■r ( 
T 


Sin 6 /O u) + ( r Sin 9 ^ u) = 0 (A. 1.20' 


u 


9u ' ^ V 
9r r 



1 9 du. 

- 5 ae ^5 ae^ 


(A. 1.21) 


and for Prandtl number equal to unity and no heat transfer 
case 

H = Hg (A. 1.22) 


Boundary conditions 

at 0=9^ u=v=0 

at 0 -*■ edge u ■* u^ 


I (A.1,23) 
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Fon-dimensionalizing the variables 


R = 


^ _ u 

IT » ^ - u 


, P = 


£. 

Pq 


(A. 1,24) 


H = 


H- 

2 % 


0 


Rei 


i 


i c. 


(e 9c) j V = 


/ Re^ 


V 

u 


~ (R^ Sin gPu) -t- 1^ (/OR Sin 0 v) = 0 


(1,1.25) 


,, 9u ^ V au 


/Or" 


9 /T du N 

'‘TT 


(A.1,26) 


H = H, 


(A.1.27) 


Boundary conditions 


at 0 = 0 u = v= 0 

at 0 •* edge u '* 1 


(A. 1,28) 


Using stream fimction 

f ^ = R^ Sin 9 p u 


t'R 


- R Sin 0pv 


(A. 1,29) 


Since Boundary layer is thin it is assumed that Sin 0 = Sin 9 
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then one gets 




) 




R 


(- 


% 


PR^ Sin 9^ pR" Sin 9 ^ PR^ Sin 9^ pR" Sin 9 


®__) 


1 

PR" 


9 

a® 


T / ?© 

j-C- — 


o R Sin 9, 


0 


(A.1.30) 


Using combined How at h similarity transformation 


H 


7} = 


and 


one gets 



/ P d® 


Tjj = 


ff + f = 0 

7)7} pm 


(A. 1.31) 


(A.1.52) 


(A.1.33) 


Boundary conditions 

at =0 f^f^O 


at 


f - 1 


(A. 1.34) 


Row Go efficient of friction 

gT gfJ 


“ t; — ^ 


o ^ 


au 


as *0=9 


=' 1 ^ '™“> 


(A. 1.35) 
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Since equations (A. 1.5) and (A. 1.34) and their corres- 
ponding boundary conditions .are same, the value of f-^^(O) are 
same for both the equations. Now comparison of equations (A. 1.19) 
and (a. 1*35) shows that value of coefficient of friction is 
same in both the cases. 

This shows that for thin similar boundary layer over a 
cone in conical inviscid flow, analysis in any of the coordinate 
system rendeis same result. 
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iPPENDES. 2 

COMPUTER PROGRAMJffi POR SOLVIHG 
DIPEERENTIAL EQUATIONS 
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